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Difference between buik materiai and deposited fiim is stiown to iiave an appreciabie inffuence 
on the Casimir force. Analysis of tlie optical data on gold films unambiguously demonstrates the 
sample dependence: the dielectric functions of the films deposited in different conditions are different 
on the level that cannot be ignored in high precision prediction of the force. It is argued that the 
precise values of the Drude parameters are crucial for accurate evaluation of the force. The dielectric 
function of perfect single crystalline gold is discussed. It is used to establish the upper limit on the 
absolute value of the force. It is demonstrated that the force between films is smaller than that 
between bulk samples mainly due to the presence of voids and electron scattering on the grain 
boundaries in the films. The minimal reduction of the force is estimated as 2% for the smallest 
separations investigated in the most precise experiments. The other sample effects can reduce the 
force further but the correction is expected to be smaller than 1%. 



I. INTRODUCTION 

The Casimir force between uncharged metaUic plates attracts considerable attention as a macroscopic manifesta- 
tion of the quantum vacuum |M IH ly| - With the development of microtechnologies that routinely allow control of 
the separation between bodies smaller than 1 fim the force became a subject of systematic experimental investigation. 
Modern high precision experiments were made using different techniques such as torsion pendulum Q , atomic force 
microscope (AFM) la Bl, rn icroelectromechanical systems (MEMS) JJ,, T3. .14] and different geometrical configura- 
tions: sphere-plate plate-plate 0|, crossed cyhnders 0. In most cases the bodies were covered with gold 
evaporated or sputter deposited to the thickness of 100-200 nm. In the only plate-plate configuration experiment 
R3| the bodies were covered with chromium. One of the bodies was covered with copper in the MEMS experiment 
|l3l . The root mean square (rms) roughness of deposited metal films was varied from 30-40 nm as in the MEMS 
experiments 0, ITsL IT^ to nearly atomically flat surfaces as in the crossed cylinders experiment [loj . Relatively low 
precision, 15%, in the force measurement was reached for the plate-plate configuration [l^l because of the parallelism 
problem. In the torsion pendulum experiment j7[ the force was measured with the accuracy of 5%. In the experiments 
E US 113 it was found with 1% precision. In the most precise up to date experiment |l3l Il4j the improvement was 
achieved due to the use of the dynamical method [llL Il2 | . Additionally the change in the resonance frequency of the 
mechanical oscillator was measured using the phase jump instead of the resonance behavior of the amplitude. In this 
way the force between gold and copper covered bodies was found with a relative accuracy of 0.25%. 

To draw any conclusion from the experiments one has to predict the force theoretically with the precision comparable 
with the experimental errors. It is a real challenge to the theory because the force is material dependent; it is very 
difhcult to fix the material properties since different uncontrolled factors are involved. In its original form, the Casimir 
force 

F,(a) = -^^ (1) 



X 

, was calculated between the ideal metals. It depends only on the fundamental constants and the distance between the 
- - 1 plates a. The force between real materials was found for the first time by Lifshitz 0,^^. The material properties 
enter in the Lifshitz formula via the dielectric function e [iC,) at imaginary frequencies lo = iC,. Correction to the 
expression is significant at the separations between bodies smaller than 1 ^m. To calculate the force the Lifshitz 
formula is used with the optical data taken from the handbooks ^3 Q| ■ The data are available only up to some 
low-frequency cutoff uJc- For good metals such as Au, Al, Cu the data can be extrapolated to the lower frequencies 
LO < LOr with the Drude dielectric function 
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which includes two parameters: the plasma frequency Wp and the relaxation frequency uJt- These parameters can be 
extracted from the optical data at the lowest accessible frequencies. The exact values of the Drude parameters are 
very important for the precise evaluation of the force. 

When two plates are separated by the distance a, one can introduce a characteristic imaginary frequency Qch = c/2a 
of electromagnetic field fluctuations in the gap. The important fluctuations are those that have frequency C ^ Cc/i- 
For a ^ 100 nm the characteristic frequency is in the near infrared region. This qualitative consideration is often 
continued to the real frequency domain fl^ 12(1 l2ll | with the real characteristic frequency uJcH = c/2a. Because 
ujch is in the near infrared it is stated that the low-frequency behavior of the dielectric function is not significant for 
the Casimir force. It is not difficult to see that this consideration is wrong [2^ I23L |2^ . The force depends on the 
dielectric function at imaginary frequencies e (j-C)- With the help of the dispersion relation e (iC) can be expressed via 
the observable function e" (w) = Ime (w): 



Because for metals e" {uj) is large at low frequencies, the main contribution in the integral in Eq. (|3Jl is given by the 
low frequencies even if ( is in the visible range. The characteristic frequency is well defined on the imaginary axis 
Cc/i = c/2a but it is wrong to introduce it in the real frequency domain. This aspect is carefully discussed in this 
paper. 

Two different procedures to get the Drude parameters were discussed in the literature. In the first approach poll25| 
the plasma frequency was fixed independently on the optical data assuming that each atom gives one conduction 
electron (for Au) with the effective mass equal to the mass of free electron. The optical data at the lowest frequencies 
were used then to find ujr with the help of Eg. In this way the parameters ojp = 9.0 eV and uJt — 0.035 eV have 
been found. In the second approach [22,|2^|23 no assumptions were made and both of the parameters were extracted 
from the low- frequency optical data by fitting them with Eq. When the best data from Ref. were used the 
result !2^ was close to that found by the first approach, but using different sources for the optical data collected in 
Ref. 15] an appreciable difference was found ||22,i22|. This difference was attributed to the defects in the metallic 
films which appear as the result of the deposition process. It was indicated that the density of the deposited films 
is typically smaller and the resistivity is larger than the corresponding values for the bulk material. However, it was 
difficult to make definite conclusions without analysis of significant amount of the data on the deposited metallic films. 
In this paper this analysis is performed and it is proven that there is a genuine sample dependence of optical properties 
of gold. First of all it means that the gold films prepared in different ways will have different Drude parameters. The 
change in the optical properties will lead to the sample dependence of the Casimir force, which is carefully investigated 
in this paper. The main ideas were outlined before |0 with the conclusion that the absolute value of the Casimir 
force between gold films is of about 2% smaller than between bulk material at the smallest separations investigated 
in the experiments [^Il4j. Here we come to the same conclusion but present much more details and facts. 

The paper is organized as follows. In Sec. ^Jit is explained qualitatively and quantitatively why the exact values 
of the Drude parameters are crucial for the precise calculation of the Casimir force. Existing optical data for gold 
are reviewed and analyzed in Sec. IIIII It is unambiguously demonstrated that there is the sample dependence of the 
optical properties of gold films. In Sec. IIVI the dielectric function of the perfect gold single crystal is discussed and 
the upper limit on the Casimir force is deduced. In Sec. a simple two parameter model is proposed allowing us 
to take into account the main reasons for the sample dependence of the Casimir force: voids and grains in the films. 
The volume fraction of voids and the mean grain size are bounded from the available optical data and the smallest 
sample correction to the force is estimated. In Sec. IVII necessarv experimental information on the films that has to 
be known for the precise evaluation of the force is briefly discussed. Comparison between the AFM experiment and 
prediction of this paper is given. Our conclusions are summarized in the last section. 
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II. IMPORTANCE OF THE PRECISE VALUES OF THE DRUDE PARAMETRS 



Let us consider the Casimir force between similar metallic plates separated by a distance a at a temperature T. It 
depends on the material optical properties via the dielectric function e (uj) at imaginary frequency ui = p^ : 



3 



Fpp (a) 



/ dqqko (r^^ exp (2afco) - l)"^ + (r^ ^ exp (2a/co) - l) "^ 



ri=0 



(4) 



where q is the wave vector along the plates [q — |q|). This formula includes the reflection coefficients for two 
polarization states s and p which are defined as 



^ _ fcp - fci ^ _ e (iCn) fcp - fci j-g-j 
fco + fci ' ^ e (iC„) fco + fci ' 

Here fco f^nd fci are the normal components of the wave vectors in vacuum and metal, respectively: 



ko^VCU^+f, fci- v/£«n)a/c2+g2. (6) 

The imaginary frequencies here are the Matsubara frequencies 



2TrkBT 

Cn = — ^ — n. (7) 

For most of the experimental configurations the force has to be calculated between sphere and plate. While the 
sphere radius R is much larger than the minimal separation a, the force can be approximately calculated with the 
help of the proximity force approximation (29J that gives the following result: 



Fsp (a) = ksTR j dqq [in (l - exp (-2afco)) + In (l - exp (-2afco))] . (8) 

This approximation is good with the precision ^ a/i? that is typically better or on the level of 0.1%. 

The material dielectric function e{iC,) is the main source of uncertainty for precise prediction of the Casimir force. 
This function cannot be measured directly but it can be expressed via the measurable function e {to) = e' (w) + ie" {uj) 
with the help of the well known dispersion relation The equation Q shows that the Casimir force is intimately 
connected with the dissipation in the material given by e" {uj) that is, indeed, the result of the fluctuation-dissipation 
theorem (see discussion in Ref. [soj'). 

The exponent exp (2afco) in Eqs. © shows that the main contribution to the force comes from the region 

a/C^ /c2 + g2 ^ i/2a. For this reason one can say that the field fluctuations with frequencies nearby the characteristic 
imaginary frequency C,ch = c/2a give the main contribution to the force. This frequency is in the infrared region for 
typical separations a = 100 — 500 nm investigated in the experiments. Therefore, for precise calculation of the force 
we have to know the dielectric function at imaginary frequencies C ' y Cr .h with the highest possible accuracy. This 
qualitative analysis is often extended to the real frequencies w 0, |23, |^ that is, of course, not true. As follows 
from the dispersion relation (|3Jl real frequencies w essentially different from c/2a dominate in e (iCch)- Suppose that 
the material can be described by the Drude dielectric function 0. Then the integrand in Eq. Q has the form 

LUpLu-riuj'^ + uJt)~^ ("-^^ + Cch) ^ ■ This fuuctiou is significant in the range uj < lut Cch- Since lut corresponds to the 
far infrared region, we have to know the function s" (w) with the highest precision in the far infrared but not in the 
infrared. Let us consider this important conclusion in detail. 

Throughout this paper the energy units will be used to measure the frequency. The following transition coefficients 
are helpful to transform it in the wavelength or in rad/s: 



1 eV ^ 1.2A1 ^im-\ 1 = 1.519 • 10^^ —. (9) 

s 

For precise calculation with Eq. Q or JSJ) the dielectric function e" (w) is usually taken from the handbooks. Fig. 
comprises the handbook data for gold The solid line shows the actual data taken from two original sources: the 
points to the right of the arrow are those by Theye [31I and to the left by Dold and Mecke [s^l . For frequencies smaller 
than the cutoff frequency ojc (0.125 eV for these data) the data are unavailable and e" (uj) has to be extrapolated into 
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FIG. 1: Handbook data for Au fv^ (solid line). The arrow separates the data from the original sources [s^l (on the left) and 
(on the right). The dotted line shows the extrapolation to low frequencies with the Drude parameters indicated in the 
upper right corner. The cutoff frequency luc and the edge of the interband transition uio separate three frequency domains 
marked as 1, 2, and 3. 



the region uj < ujc- This is a dehcate procedure because the contribution of the low frequencies dominates in e (iC)- 
In the field of the Casimir effect the following procedure has been adopted (25| . The plasma frequency is calculated 
using the relation 



= 7^ (10) 

where N is the number of conduction electrons per unit volume, e is the charge and to* is the effective mass of 
electron. Assuming that each atom gives one conduction electron, that the effective mass coincides with the mass of 
free electron and using the bulk density of Au it was found LOp = 9.0 eV . With this value of the plasma frequency the 
best fit of available low-frequency data with the Drude function (0) gave for the relaxation frequency lOt = 0.035 eV 
[25l |. The dotted fine in Fig. ^ shows the Drude extrapolation with these parameters. 

One can separate three frequency regions in Fig. ^marked as 1, 2, and 3. In the high energy domain uj > loq, where 
ujQ ^ 2.45 eV 31] is the edge of interband absorption, the Casimir force only weakly depends on the behavior of the 
dielectric function. However, as we will see in the next section, analysis of the interband absorption can give valuable 
information about voids in the metallic films. The region 2 extends from the cutoff frequency ujc to the absorption 
edge loq. The curve in this region is also sensitive to the sample quality. It gives perceptible contribution to e (iC) but 
what is more important is that this domain completely defines the Drude parameters. These parameters will be used 
to extrapolate the curve to the low-frequency region 1, uj < Wc- The low-frequency region gives the main contribution 
to £ (iC) and for this reason the precise values of the Drude parameters are very important. 

Let us write the dielectric function at imaginary frequencies in the following form: 



e «) = 1 -I- £i (zC) + 62 {iO + 63 (iC) , (11) 

where the indexes 1,2 and 3 indicate the integration range in Eq. ||2Jl that are < ^ < Wc, Wc < < woj Siiid 
ujo < ijJ < oo, respectively, ei can be calculated explicitly using the Drude function (01. It gives 



£i «) 



tan 



c 



■ tan 



(12) 



Two other functions 62 and £3 have to be calculated numerically. The result for all functions including e (i^) is 
shown in Fig. [3 One can clearly see that £1 {iQ dominates in the dielectric function at imaginary frequencies up to 
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FIG. 2: Contributions of different real frequency domains to the dielectric function on the imaginary axis e{i(). The low- 
frequency {uj < u)c) contribution ei(jC) clearly dominates while C < 5 eV. 



C ~ 5 eV. It means that for all experimentally investigated situations the region 1, where the extrapolated optical 
data are used, gives the main contribution to e(i(^). For example, for the separation between bodies a — 100 nm 
the characteristic imaginary frequency is Cch — 0-988 eV. At this frequency the contributions of different domains to 
£ (iCch) are ei — 68.42, 62 = 15.65, and £3 = 5.45. 

In the following section it will be demonstrated that the dielectric function of the films used for the force mea- 
surement can deviate considerably from the handbook data. The main reason for this is genuine sample dependence 
of the optical properties. Here let us consider the question how sensitive is the Casimir force to this deviation. To 
this end consider small variations Ae„ of the dielectric function £„ = e(i^„). The relative change in the force can be 
written via the relative variation of e„ as 



(13) 



where F„ is the n-th term in the sum Q or (jSJ and F is the total force. To see how important is the change of the 
dielectric function at real frequencies, it is more convenient to introduce instead of Asn/sn the relative variations for 
each frequency domain i, ^, or 3: 

^ = ^A„r + ^A„. + ^A„3 (14) 

with 



A„j — 



7re 



nj JD 



ujAe" {lj) 



J = 1,2,3. 



(15) 



Here Snj is defined as Sj (iCn), Ae" {uj) is the measured deviation at real frequencies and Dj is a particular frequency 
interval coinciding with the domain j. With these definitions the relative change in the force can be presented as a 
linear combination of A„j : 



AF 



^ 00 



n=0 



\ dSn 



Snl A„i 



dFa 

dSn 



A, 



dFn 
de. 



in3 



(16) 



The coefficients in front of Aj can be called the weight functions. They describe the weight with which domain j 
contributes to the force. These weight functions are shown in Fig. |21for the separation a = 100 nm. One can clearly 
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see that the weight for the domain 1 is the largest. When a increases the weight functions become more narrow and 
the domains 2 and 3 become less important. In the next section the relative variation of e„j for real gold films will 
be discussed where more detailed information on A„j will be presented. 

The zero term n = in Eqs. I|13|) - (|16ll h as to be described separately. In the Lifshitz formula or (|SJ) the n = 
term rose controversy in the literature [22, l33L IsiL ISa . 1371 133. l39l W\ . It is still an unresolved problem in the field of 
the Casimir effect. The point of disagreement between different authors is the value of the refiection coefficient in 
the limit C ^ 0- Bostrom and Sernelius found rg = applying the Drude behavior to e {iQ at C — > 0. Svetoyoy 
and Lokhanin "25, "27, "371 got = 1 demanding continuous transition to the ideal metal at C 0. Bordag et al. |34| 
have used the dielectric function of the plasma model in the limit C ^ ^iiid found 




There is no problem with Vp coefficient which is always rp = 1 in the C ~* limit. 

In the first two approaches and Vp do not depend on the material and the zero term should be omitted in Eqs. 
(|13ll - (|16|l . In the plasma model approach Tg depends on the plasma frequency and the n — term will change 
when ujp varies. The contribution of the zero term to Eq. will be 



\f). era term [du;p"'~' 



u^p ] (18) 



In Fig. 13 the coefficient in front of Aujp/ujp is shown by the dot on the ordinate axis. 

The main conclusion that can be drawn from the analysis above is that the deviation of the Drude parameters ujp 
and uJt from their handbook values gives the largest contribution to the change of the Casimir force. 



III. OPTICAL DATA FOR GOLD AND ANALYSIS 



In this section we are going to show that the optical properties of gold films depend on the details of the deposition 
method and following treatment to the extent that cannot be ignored in the calculation of the Casimir force. Optical 
properties of gold were extensively investigated in 50- 70th. Today most of the researches are interested in gold 
nanoclusters but these data are inappropriate for our purposes. There are only a few exceptions. Reflectivity of 
evaporated gold films was investigated |41| in the wavelength range 0.3 — 50 /iw to analyze its dependence on the 
mean grain size. Since only normal reflectance has been measured, one cannot get the information on both real and 
imaginary part of the dielectric function but these data are very helpful. The reflectivity of gold in submillimeter 
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FIG. 4: The imaginary part of the dielectric function in the interband region for different samples. The sohd hne represents 
well annealed bulk-like film by Theye [sj. The dots are the data by Johnson and Christy 14^ found for unannealed films. The 
dashed and dash-dotted lines are the modern data by Wang et al. |43 | for unannealed films. They correspond to the films 
deposited with e-beam and thermal evaporation methods, respectively. 

range, 513 /im, was measured with high precision f42i |. In Ref. [43j | the films evaporated on glass with two different 
methods were characterized ellipsometrically in the photon energy range 1.5 — 4.5 eV. This work allows one to compare 
the interband absorption with the older measurements. In addition, the optical properties of films deposited with the 
magnetron dc sputtering and with the filtered arc process ^44] were measured at the photon energies 1.5 — 3.5 eV and 
compared with each other. One can indicate also two recent studies |45l where new ellipsometric equipment was 
tested but the preparation of the films was poorly described. We will use the recent papers where it is possible but 
most of the data will be taken from the old works. In many old works the importance of sample preparation methods 
was recognized and carefully discussed. Complete bibliography of the old publications up to 1981 one can find in Ref. 

A. Interband spectral range 

Let us start from the interband absorption or domain 3. This range is not very significant for the Casimir force 
but it provides us with important information how the sample preparation procedure influences the data. There is 
significant amount of data in this range obtained by combined reflectance and transmittance l3ll [4 p| , e llipsometric 
spectroscopies "i^, H on unannealed or annealed ElEiii thin film [siTlS l itl lil IST f 

or bulk 49, 50] samples measured in air 43, 44, 48, 57] or ultrahigh vacuum |3lL llgllSOl ] . A representative picture of 
typical variations of e" {oj) is shown in Fig. ^ The data for this figure were chosen to represent very dense (bulk-like) 
film |31|], the films similar to those used for the Casimir force experiments and some intermediate case [i^ . 

Theye [sj] described her films very carefully. The samples were semitransparent Au films with the thickness 
100 — 250 A evaporated in ultrahigh vacuum on supersmooth fused silica. The deposition rate was a few A/s and the 
substrate was kept in most cases at room temperature. After the deposition the films were annealed under the same 
vacuum at 100— 150° C, the process being monitored by simultaneous dc resistance measurements. The structure of 
the films was investigated by X-ray and electron-microscopy methods. The "good" films were continuous formed of 
joined regularly shaped crystallites, the average lateral size of which ranged from 300 to 500 nm. The main defects 
were grain boundaries. It was demonstrated that flatness and smoothness of the film surfaces were very good. The 
dc resistivity of the films were found to be very sensitive to the conditions of preparation. The errors in the optical 
characteristics of the films were estimated on the level of a few percents. 

Johnson and Christy |48| evaporated their films onto fused quartz substrates at room temperature. The evaporation 
rate was relatively high, 60 A/s, and the thickness was in the range 250-500 A. The presented data were found 
for unannealed samples. The structural characterization of the films was not made. Errors in the optical data were 
estimated on the level of 2%. 

Wang et al. 113] evaporated the films on fused silica or optical glasses with two different methods. For the thermal 
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evaporation the rate was 50-100 A/s and for the electron-beam evaporation it was 5 A/s. The substrates were kept at 
room temperature and there was no annealing after the deposition. Roughness of the films or structural details were 
not reported. Opaque 150 nm thick films were characterized ellipsometrically in the photon energy range 1.5 — 4.5 eV. 
Errors were not specified but typical precision of the method is of about 1%. 

As one can see in Fig. 0]the spectra differ mainly on the scale factor. The experimental errors are too small to explain 
significant variation of the data. One can conclude that there is genuine sample dependance of optical characteristics 
of gold films. Differences among these spectra have been attributed to various effects such as roughness Isil Isi l . surface 
films 52], unspecified internal electron-scattering defects _31, 53], strain [HJ], grain-boundary material |5ll Is^. IssI]. 
and voids ,56] . Special investigation directed to quantitative understanding the discrepancies on a systematic basis 
was undertaken by Aspnes et al. [s^] (see also [s^] and H^)- The authors analyzed the literature data and their 
own gold films evaporated by electron beam or sputter-deposited. The best films were obtained on cleaved NaCl 
substrates. Glass, Si02, and sapphire substrates were also investigated. Air cleaved NaCl substrates were superior 
because gold start to grow on them epitaxially. Scaling behavior of the spectra was attributed to different volumes of 
voids in the films prepared by different methods. It was also stressed that grains in a film imply voids, because it is 
impossible to produce grain boundaries in fee or hexagonal close-packed lattices without losing density. 

To describe the contribution of voids in the dielectric function, the Bruggeman effective medium approximation 
|60l] was used |53]. Nowadays this approximation is successfully applied for interpretation of the spectra of different 
nanocomposite materials (see, for example, |6ll Ig^V If e is the dielectric function of Au in its homogeneous form, 
then the effective dielectric function, (e), of the material containing a volume fraction of voids fy can be found from 
the equation 



where Eh is the dielectric function of the "host" material. In the Bruggman approximation it is assumed that the 
host material coincides with the effective medium, en = (e), so it treats both void and material phases on an equal. 
If the void fraction is small, fy <^1, then the solution is 



(e>-.(l-3/.^). (20) 

Aspnes et al. jsj analyzed available data for e" {uj) in the interband spectral range fitting them with Eq. H19|) . It 
was found out that in their e-beam evaporated unannealed sample the volume fraction of voids was 4% larger than 
in Theye annealed thin films [sj but 5% smaller than /y in Johnson and Christy |_42j unannealed thin films. One of 
the films was sputter-deposited on NaCl substrate that was held at liquid nitrogen temperature. This sample was 
full of voids that were clearly seen with the transmission electron microscopy. The volume fraction was estimated as 
fv = 25% in respect to Theye 31j film. One should mention that roughness of the film also will reduce e" [lo) in the 
interband region but this is a minor effect (53 because the films in Refs. |^ |53 were quite smooth. However, 
significant roughness could be the reason why Wang et al. thermally evaporated film shows low absorptivity. 

Concluding this discussion one can say that the sample preparation procedure significantly infiuences the optical 
properties of the films in the interband spectral range. At first sight it seems unimportant for the Casimir force since 
variation of the dielectric function in this range gives only minor correction to the force. However, the volume fraction 
of voids, fv, that can be estimated from these data defines the exact value of the effective plasma frequency which is 
very important for the precise prediction of the Casimir force. To see it clearly consider Eq. (I20II at low frequencies 
where ]e] ^ 1. In this limit the effective dielectric function will be simply proportional to that of the homogeneous 
material 

(21) 

In the Drude region where e is given by Eq. 121 this relation is equivalent to the existence of the effective plasma 
frequency 

{u.l)^ul{l-\fy^ (22) 
because voids reduce the number of conduction electrons per unit volume. 
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Correction to the force due to variation of e" (w) in the interband spectral range is expected to be small. Neverthe- 
less, let us estimate it. For numerical calculation it was assumed that the bulk material can be described by Theye 
data but the film material is given by Johnson and Christy data Besides, it was assumed that the difference 
between these data is negligible at the photon energies larger than 6.5 eV. For the third term in Eq. Hlt)|l and for the 
sphere-plate geometry it was found (AF/F)^ = —0.004, —0.002, —0.001 for the separations a — 50, 100, 150 nm, 
respectively. Indeed, the correction is negative since any deviation from the perfect material will reduce the reflection 
coefficients and, therefore, the absolute value of the force will be smaller. At separations larger than 150 nm the 
correction to the force is completely negligible on the level of 0.1%. 



B. Intraband spectral range 



When the photon energy is below the interband threshold, uj < loq, only free carriers should contribute to the 
dielectric response. However, the Drude function ^ obviously fails to describe the data nearby the transition (see 
points below lvq in Fig. 0)). Steep rise of e" near uq = 2.45 eV originates from d bands to Fermi-surface transition near 
L point but there is a long tail which extends well below 2 eV. This tail cannot be explained by the temperature 
broadening of the threshold. It was demonstrated [H^l that an additional transition at 1.94 eV near X point in the 
Brillouin zone is involved. This transition has a smooth edge due to topological reason. In addition, the surface 
roughness and the grains can contribute to the tail broadening [H^ . This tail is not very important for the Casimir 
force but it marks the region where the Drude parameters cannot be determined precisely. 

In contrast with the interband region there are only a few sources where the dielectric function was measured in 
the infrared. The available data for e" (oj) in the range uj < loq are presented in Fig. The dots represent the points 
from the handbook 17] comprising Dold and Mecke data [s^l for w < 1 eV and Theye data l3j| for higher frequencies. 
Dold and Mecke did not describe carefully the sample preparation. It was reported only that the films were 
evaporated onto a polished glass substrate and measured in air by using an ellipsometric technique. Annealing of the 
samples was not reported. These data are typically referred in connection with the Casimir force calculation. The 
squares are the data by Weaver et al. jS^] that were found for the electropolished bulk Au{110) sample. Originally 
the reflectance was measured in a broad interval 0.1 < a; < 30 eV and then the dielectric function was determined 
by the Kramers-Kronig analysis. Due to indirect determination of e the recommended accuracy of these data is only 
10%. Motulevich and Shubin result for Au films is marked with the open circles. In this paper the films were 
carefully described. Gold was evaporated on polished glass at pressure ~ 10~^ torr. The investigated films were 
0.5 — 1 fj-m thick. The samples were annealed in the same vacuum at 400° C fore more than 3 hours. The optical 
constants n and k {n + ik = y/e) were measured by the polarization methods in the spectral range 1 — 12 /im. The 
errors in n and k were estimated as 2-3% and 0.5-1%, respectively. Moreover, the density and conductivity of the 
films were measured. The crosses represent Padalka and Shklarevskii data [gsI ] for evaporated unannealed Au films. 
The solid line shows the only modern data for Au films sputtered on Si substrate. No additional details about 
the film preparation were reported. The precision of optical data was claimed to be better than 1% in the entire 
spectral range. 

Different researches stressed that the conduction electrons are much more sensitive to slight changes in the material 
structure [sil IstI l6^ . As one can see in Fig. [S] it is, indeed, the case. Variety of the data in any case cannot 
be explained by the experimental errors. The observed scattering of the data is the result of different preparation 
procedures and reflects genuine difference between samples. Bennett and Ashley |H3 observed considerable reduction 
in the reflectance of Au films when instead of ultrahigh vacuum (5 • 10~^ torr) the films were evaporated in the ordinary 
vacuum (10^"' torr). It was noted also that it was more difficult to prepare films with reproducible reflectance in the 
ordinary vacuum. Even more important are the deposition method (there is a whole spectrum of evaporation and 
sputtering methods), type of the substrate, its temperature and quality, and deposition rate. When we are speaking 
about the Casimir force measurement with the precision of 1% or better there is no any more such a material as gold 
in general. There is only a material prepared in definite conditions. 

The data in Fig. |31 below ojq should be described quite well with the Drude function 



e"{u;) = 



UJ ((jj2 -f 0^2) 



(23) 



While w 3> Wr the data on the log- log plot should look like a straight line with the slope —3 shifted along the ordinate 
axis mostly due to variation of ujr for different samples. The data in Fig. [S] are in general agreement with the 
expectations excluding the solid line. The latter cannot be described with the Drude behavior, probably, due to poor 
preparation of the film. Padalka and Shklarevskii data (crosses) demonstrate a long Drude tail that can be observed 
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FIG. 5; Available optical data in the infrared region. The dots represent the Dold and Mecke data at oj < 1 eV for evaporated 
unannealed films ^5] and Theye data [sH for higher frequencies. The squares are the data [4^ for electropolished bulk ^m(IIO). 
The open circles are the data for carefully prepared well annealed evaporated films |6^ . The crosses represent the data ;65.] for 
unannealed evaporated films. The only modern data 46] are given by the solid line. 



up to 0.3 eV. To all appearance it is again the result of poor sample preparation. The Drude parameters can be 
found fitting both e' and e" with the function ((2Jl. This procedure will be discussed later. 

Now let us estimate the change in the Casimir force due to variation of e" (uj) in the infrared. For this purpose 
we calculate A„2 in Eq. H15|l using as Ae" (w) the difference between the squares 47| and points in Fig. [S] in 
the interval ojc < uj < ujq. With this A„2 the second term in Eq. (|16|l is calculated that gives the relative change 
in the Casimir force due to variation of s" {lu) in the indicated frequency interval. For the sphere-plate geometry it 
was found that (AF/F)^ = -0.025, -0.020, -0.016 for the separations a = 50, 100, 150 nm, respectively. The 
change in the force is negative and quite large. The main reason for this is the difference between the relaxation 
frequencies of the bulk material (squares) and poorly evaporated film (dots). As one can see in Fig. El below 1 eV 
the dielectric function s" for the film is roughly two times larger than that for the bulk material (squares) so A„2 is 
approximately —0.5 for all important n. There is no contradiction with the expectation that the force must be larger 
for the bulk material. We took into account here only the frequency interval < ix) < cjq- When the contribution of 
the low-frequency domain < ui < uic will be added the total force will be larger for bulk gold. 



C. The Drude parameters 



The dielectric function in the low- frequency domain, lu < ujc, where direct optical data are inaccessible, is defined 
by the dielectric function in the infrared domain, lUc < uj < ujq, because one has to extrapolate the optical data from 
this domain to inaccessible region. The real and imaginary parts of e follow from the Drude function Q with an 
additional term 7^ in e'; 



s'H = r-^^, s"H^ /f7 (24) 

The polarization term V has appeared here due to the following reason. The total dielectric function e = e(^c) + ^(i) 
includes contributions due to conduction electrons e(c) and the intcrband transitions £(.;). The polarization term 
consists of the atomic polarizability and polarization due to the intcrband transitions e'^^^ 



^ = l + ^+4oH' (25) 

where a is the atomic polarizability and Na is the concentration of atoms. If there is no the Drude tail the last term 
can be considered as a constant. This is because the intcrband transitions have the threshold behavior with the onset 
ujQ and the Kramers-Kronig relation allows one to express e',^^ as 
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N 


ujp, (eV) 


ujt ■ 10^ (eV) 


r 


Ref., mark 


1 


7.52 ±0.21 


6.1 ±0.9 


-28 ± 67 


[17,32], ■ 


2 


8.41 ±0.01 ±0.42 


2.0 ±0.04 ±0.3 


7.3 ±6.0 


mj' ■ 


3 


8.79 ±0.12 ±0.09 


4.2 ±0.3 ±0.2 


14 ±79 


[64], O 


4 


6.81 ±0.07 


3.6 ±0.3 


-21 ±39 


[65], X 


5 


9.04 ± 0.03 


2.67 ±0.03 


24 ± 117 


[68] 



TABLE I: The Drude parameters found by fitting the available infrared data for e' (lj) and e" (uj) with Eq. 124^ . The first 
error is the statistic one and the second error (if present) is the systematic one. 



= - dx^^. (26) 

For LU ^ luq this integral does not depend on uj. In reality the situation is more complicated because the transition is 
not sharp and many factors can influence the structure of the Drude tail. It will be assumed here that V is a constant 
but the fitting procedure will be shifted to frequencies where the Drude tail is not important. In practice Eq. H24|) 
can be applied for uj < 1 eV. 

As we have seen above the plasma frequency ujp of a sample under investigation can deviate from that of perfect 
bulk material due to presence of voids (see Eq. (|22|l ). The Casimir force is measured between films which are far 
from perfect. The best optical properties demonstrate well annealed films prepared in ultrahigh vacuum |3lLl66l |. The 
annealed films were not used for the force measurement and it is hardly reasonable to do in any future experiments 
|90j | . For high precision calculation of the force the procedure proposed in Ref. 25] is not good any more. One cannot 
use the bulk plasma frequency to extrapolate the dielectric function to experimentally inaccessible frequencies. The 
second Drude parameter, the relaxation frequency tOr, is even more sensitive to the sample preparation procedure. 
Therefore, the Drude parameters have to be extracted from the optical characteristics of the same sample which is 
used for the force measurement. Our purpose here is to establish the magnitude of the force change due to reasonable 
variation of the optical properties. To this end the available low-frequency data for e" (ui) shown in Fig. [Sjand e' (cu) 
(not shown) were fitted with Eq. (|24|) . The results together with the expected errors are collected in Table UJ 

The first error in Tabled is the statistical one. It was found using criterion for joint estimation of 3 parameters 
|69j |. The second error is the systematic one. It was found propagating the optical errors, when indicated, to uip and uJr- 
One can see that the polarization term V cannot be resolved in the frequency range where e' is very large. Significant 
variation of the plasma frequency, well above the errors, is a distinctive feature of the table. Carefully prepared 
annealed samples demonstrate larger value of ujp. The rows 1 and 4 corresponding to the evaporated unannealed films 
have considerably smaller ujp. Note that our calculations are in agreement with that given by the authors [s^ 
themselves. To all appearance these low values of ujp are the result of poor preparation of the films. In accordance 
with the effective medium approximation at low frequencies the effective plasma frequency is given by Eq. I|22|) . Using 
this equation and uip — 9.0 eV for homogeneous gold the fraction of voids fv — 0.20 and 0.28 was found for the rows 
1 and 4, respectively. Discussing the interband transitions we have seen that this large fraction of voids is possible 
for poor prepared films like that sputtered at liquid nitrogen temperature [s^ or thermally evaporated ^43,] . 

The fifth row in Table ^should be discussed separately. It is based on the reflection data [63 in the infrared for 
carefully prepared opaque films evaporated and annealed in ultrahigh vacuum on extremely smooth fuzed quartz 
substrate. The films have the highest reflectivity in the infrared of that ever reported. These data have been analyzed 
by Bennett and Bennett |68l | using the conductivity of bulk material and assuming that there is one conduction 
electron per atom. Perfect agreement with the Drude theory in the wavelength range 3 — 32 fim was found with no 
adjustable parameters. The errors in the row 5 demonstrate rather the level of this agreement than the statistical or 
systematic errors. 

To have an idea how good is the fitting procedure, in Fig. the experimental points and the best fitting curves 
are shown for Dold and Mecke data (dots and solid lines) and Motulevich and Shubin data [g^I (open circles and 
dashed lines). One can see that for e" at high frequencies the dots lie above the solid line indicating the presence of 
the Drude tail. Coincidence of the solid and dashed lines for e" is accidental. The fits for e' are nearly perfect for 
both sets of the data. 

Calculating the Casimir force Lambrecht and Reynaud ^25*] assumed that the plasma frequency can be taken for 
bulk Au with one conduction electron per atom and found ujp = 9.0 eV. In contrast with '6Sj this assumption was not 
supported by any optical data. The relaxation frequency ujr = 0.035 eV was found then by fitting Dold and Mecke 
data [33 for e". The curves corresponding to this set of parameters are shown in Fig. Elas dotted lines. It is clear 
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FIG. 6: The infrared optical data by Dold and Mecke (dots) and by Motulevich and Shubin |64|l (open circles) together 
with the best Drude fits given by the solid and dashed lines, respectively. The dotted lines present the fit which is typically used 
in the calculation of the Casimir force |25ll . The latter one much better agrees with the open circles than with the handbook 
data (dots). 



that there is a significant mismatch for e'. The dotted lines much better agree with Motulevich and Shubin data |6J| 
than with the data they should describe. If the force is calculated with the handbook data, the extrapolation to low 
frequencies has to be done with the parameters which are in agreement with the optical data. These parameters are 
given in the row 1 (Table P). 

Let us estimate how sensitive is the force to the variation of the Drude parameters. To this end the first term in Eq. 
()16|l is estimated. In the domain of interest, < uj < lUc, we calculate Ae" (w) in Eq. H15|l as the difference between the 
function H23|) with ujp — 7.52 eV, uJt — 0.061 eV (evaporated film, first row in TableP) and the same function with the 
bulk-like parameters Wp = 9.0 eV, u)r = 0.035 eyjU. The result is the following: {AF/F)^ = -0.070, -0.063, -0.054 
for the separations a = 50, 100, 150 nm, respectively. The correction is negative since the evaporated film has smaller 
reflectivity in comparison with the bulk material. Its magnitude is large because the force is more sensitive to the 
behavior of e" (lu) in the far infrared than in the infrared or visible range as was discussed in Sec. ^ At any 
circumstances this large correction cannot be ignored. 

The main problem with precise calculation of the Casimir force is to find the correct values of the Drude parameters 
to extrapolate the dielectric function to low frequencies. Obviously, using the bulk material parameters we overestimate 
the force which is measured between the deposited films. Discussion in this section demonstrates that the film quality 
is improved after annealing. Unannealed films used for the force measurement will inevitably contain some fraction of 
voids independently of the used equipment. This fraction can be larger than 10% as for the thermally evaporated films 
|43| but, probably, not smaller than 4% as for the films evaporated with e-beam on NaCl substrate [5^ or sputtered at 
room temperature [4^ Is^ . The void fraction will influence the effective plasma frequency. The relaxation frequency 
is even more sensitive to the film quality because of electron scattering on the grain boundaries and other defects. 
To be able to predict the force with the precision better than 1% the films have to be carefully characterized. 

At the end of this section let us estimate the contribution of the Drude tail to the force. It is calculated in the 
frequency range 1 eV < uj < loq that was not used in the determination of the Drude parameters. We calculate 
Ae" (uj) in Eq. (|15() as the difference between e" given by Eq. (|23|l with the parameters in the row 1 (Tabled and 
e" from the handbook data. For the contribution to the force between sphere and plate it was found {AF/F)^ — 
—0.003, —0.002, —0.001 for the separations a = 50, 100, 150 nm, respectively. The effect is minor. 



IV. UPPER LIMIT ON THE CASIMIR FORCE 



While the Drude parameters of the films are not fixed precisely from the optical characterization, one can establish 
a reliable upper limit on the force. It is based on the observation H^H^I that the reflection coefficients |rs^p| in the 
Lifshitz formula (0J or © are always larger for perfect material than for any real material containing different kind 
of defects. As the result the force between perfect bodies will be the largest. To be more precise, at a fixed imaginary 
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frequency \rs\ and \rp\ are monotone increasing functions of e (zC). The dielectric function e (zC) can be represented 
in the following form: 



oo UJQ 



where e'^ (w) is the Drude function (|23|l . The first two terms in Eq. (|27() will be maximal for a single-crystal when 
ujp is the largest, cOr is the smallest, and the interband {to > loq) absorption is maximal. The same is not true in the 
range of the Drude tail, uji < oj < uiq. As one can see in Fig. |51the smallest Drude tail shows the single crystalline 
sample, while the films demonstrate much larger tail. However, the increase of e (iC) due to the tail is always small in 
comparison with the decrease of e (i^) due to deviation of the Drude parameters from the bulk values. Hence, e (i^) 
will be the largest for single crystalline gold. 

It would be natural to take the dielectric function of the single crystalline sample to find the upper limit on the 
Casimir force. This possibility was rejected due to the following reason. Precision of these data is only 10% because 
the Kramers-Kronig analysis was used to get both of the optical constants. This is the reason for large systematic 
errors in the Drude parameters (see row 2 in Table P). Since the force is very sensitive to the value of uJp, the precision 
of the data is not good for our purpose. To get the Drude parameters of the perfect sample, the plasma frequency was 
calculated with the help of Eq. (|10|l assuming that each atom gives one conduction electron with the effective mass 
m* = nie- In this way it was found Up = 9.03 eV. The relaxation frequency was calculated then from the resistivity 
of bulk gold at room temperature p = uJr/£o^p = 2.25 /ifi • cm that gave ujr = 0.025 eV. The deviation from the 
parameters presented in the row 2 is a little bit larger than one standard deviation but acceptable. Note that the 
parameters found in this way agree very well with Bennett and Bennett data (row 5 in Table H)) . 

The upper limit on the force is calculated using the Drude function H23|l at w < wi = 1 eV with the parameters 
UJp — 9.03 eV, uJr = 0.025 eV. Starting from uj = ujo = 2.45 eV Theye data 0, |^ are used as the best data 
known in the interband region. At frequencies uj > 6 eV the handbook data are used. The force is not sensitive 
to the reasonable variation of the data due to the sample effect at w > 6 eV. The most problematic interval is the 
range of the Drude tail, uji < lu < ujq. In this interval the data are interpolated linearly on log-log plot between 
e" {ui) and e" {llIq). Of course, any real data should lie below the interpolated line. The only justification for this 
procedure is the fact that we are trying to find the upper limit on the force. This interpolation gives the largest 
contribution. If, for example, the Drude function is continued up to = luq then the relative change in the force 
will be AF/F = -0.005, -0.003, -0.002 for a = 50, 100, 150 nm, respectively The function e" (uj) used for the 
calculation of the upper limit is shown in Fig. [7| by the solid line. Note that it is not just an abstract line for the 
best gold sample. In the interval 0.039 < uj < 0.41 eV it is very close to the data for the films carefully prepared in 
ultrahigh vacuum [g^- The handbook data presented in the same figure by the dashed line for u) > 0.125 eV. Two 
extrapolations to low frequencies are also shown. The dashed line at w < 0.125 eV corresponds to the best fit of 
the data in the infrared range 0.125 < uj < 1 eV providing smooth continuation of both e' (w) and e" (uj) to lower 
frequencies. The dotted line gives the extrapolation used in Ref. with the parameters ujp = 9.0 eV, uir = 0.035 eV. 
This extrapolation is continuous only for e" (uj). Mismatch of e' at uj ^ 0.125 eV is 100% as one can see in Fig. 
(dotted line). 

One can ask why in the infrared the solid line representing the best sample lies below the handbook data. Larger 
e" seems augments e (iC) and, therefore, increases the force. In reality the integral effect is important. Larger values 
of e" in the infrared happen mostly because of larger uJt (see Eq. ((221 )■ In the far infrared the function with larger 
uJt will be smaller than that for the best sample. Because the contribution of low frequencies is more important in 
the integral £ (iC) always will be larger for the perfect single-crystal. 

The results of calculation of the upper limit on the Casimir force are presented in Table ^ It gives the so-called 
reduction factors [2^ which are defined as 

_ Fpp (a) _ TT^ he 

"^""-F^y ^^^'("^-"240^ 



for the plate-plate geometry and 



_ (a) (29) 
F^piaY ^^P^"'- 360 a3 ^ 

for the sphere of radius R and plate. Here Fpp and Fgp are the forces given by the Lifshitz formulas l@J and ||SJ|, 
respectively, and Fpp, F^p are the corresponding expressions for the Casimir force between bodies made of the ideal 
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FIG. 7: The dielectric function used to calculate the upper limit on the Casimir force (solid line). The dashed line represents 
the handbook data fl7)l extrapolated to low frequencies oj < 0.125 eV continuously for both e' and e". The dotted line shows 
the extrapolation used in Ref. I25II which is continuous only for e". 



a (nm) 


Vpp 






Tjsp 




^nlp 


60 


0.364 


-0.002 


0.007 


0.435 


-0.005 


0.009 


80 


0.424 


-0.003 


0.009 


0.497 


-0.008 


0.011 


100 


0.472 


-0.005 


0.010 


0.546 


-0.011 


0.012 


120 


0.513 


-0.007 


0.011 


0.586 


-0.015 


0.013 


140 


0.547 


-0.010 


0.012 


0.620 


-0.019 


0.013 


160 


0.578 


-0.012 


0.013 


0.648 


-0.023 


0.014 


180 


0.604 


-0.014 


0.013 


0.673 


-0.027 


0.015 


200 


0.627 


-0.017 


0.014 


0.694 


-0.031 


0.015 


220 


0.648 


-0.020 


0.014 


0.713 


-0.035 


0.015 


240 


0.667 


-0.022 


0.014 


0.730 


-0.039 


0.016 


260 


0.684 


-0.025 


0.015 


0.745 


-0.044 


0.016 


280 


0.699 


-0.028 


0.015 


0.759 


-0.048 


0.016 


300 


0.713 


-0.030 


0.015 


0.771 


-0.052 


0.016 



TABLE 11: The reduction factors for the plate-plate geometry ripp and sphere-plate geometry rjsp corresponding to the upper 
limit on the Casimir force. /S.r]^^ gives the temperature correction according to Ref. |33|| and Ar;^^ gives the same correction 
according to Ref. |3^ . 

metal. The reduction factors are given for the case when the temperature correction is negligible. It corresponds to 
the plasma model prescription (34] for the n = term in the Lifshitz formula. Corrections to the reduction factors 
are also shown for the cases when the zero term is described by the Drude model Arj^^ or coincides with the 
ideal metal Arj^^ [s^- They give the temperature corrections in different approaches. 

The values presented in Table ^] one can compare with the result by Lambrecht and Reynaud |25j who gave for 
a — 100 nm the values r]pp = 0.48 and rjsp = 0.55. They are very close to that presented in Table ITTl This coincidence 
is because the increase of the force due to the infrared region was compensated by its decrease in the far infrared 
since larger value uJr = 0.035 eV was used. Therefore, the result of Ref. [2^ represents rather the upper limit on 
the force than the real force corresponding to the handbook data. The latter force should be calculated with the 
Drude parameters given in Tabled (row 1). For the sphere-plate geometry at a = 100 nm it is 6.3% smaller as 
was demonstrated in Sec. IIII CI For small separations presented in Table Hll At]/ rj gives the relative temperature 
correction. It is well known that this correction is negative for the Bostrom and Sernelius approach "3^ making 
the absolute value of the force smaller; it is positive when the n = term is assumed to coincide with that for the 
ideal metal [s^- In the latter case the relative correction is nearly constant at small separations and is of about 2% 
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independently on the geometry. In the former case it decreases from -0.5% to -4.2% for the plate-plate and from -1.1% 
to -6.7% for sphere-plate geometries. 

We did not include in the upper limit the correction due to surface roughness [iS 1^ for obvious reason: it should 
be specified for a definite experiment. This correction always increases the absolute value of the force. Additionally, 
the nonlocal effects in the interaction of electromagnetic field with a solid were not included in the upper limit on the 
force. It was demonstrated [z^l that the nonlocal screening effect due to excitation of plasma oscillations in metals 
will reduce the Casimir force on a few percents at small separations. The influence of the anomalous skin effect on 
the force was also analyzed recently jTSj. It was found that this effect gives only minor correction which reduces the 
force smaller that 0.5%. 



V. EXPECTED MAGNITUDE OF THE SAMPLE EFFECT 



As we have seen above the main uncertainty of the Casimir force comes from the uncertainty of the Drude param- 
eters. The reason for variation of the effective plasma frequency, Wp, is the presence of voids in the deposited films as 
was discussed in Sec. IIII Al The relaxation frequency, Wr, is much more sensitive to the details of the film preparation 
[31I Is^ Ig^ . It is well known that the resistivity p of deposited films is larger than that for the single crystalline 
material. The resistivity of unannealed medium thick (150— 200 nm) gold films was ranged from 4 to 10 /ifi • cm 
in dependence on the method (electron beam, sputtering) and rate (1 A/s, 20 A/s) of deposition on S'i-substrate. 
For well crystallized annealed film deposited in ultrahigh vacuum Theye |31| has found p = 2.97 /ifi • cm but for 
unannealed film this value was as high as 4.40 /ifi • cm. These values should be compared with the resistivity of bulk 
gold 2.25 p^ ■ cm. The difference between film and bulk values was attributed to the presence of grain boundaries 
in the films (sj- In this section a simple two parameter model will be presented allowing one to estimate the force 
between films via the volume fraction of voids /y and the mean grain size D. 



A. Optical conductivity as a function of the grain size 

Mayadas and Shatzkes [t^ proposed the first model for calculation of the dc conductivity, a = 1/p, of films taking 
into account the electron scattering on the grain boundaries. This model successfully describes films and nanowires 
In the model the films were treated as joined grains with bulk material properties; the grain boundaries 
contributed to the electron transport. A semiclassical approach based on the Boltzmann transport equation was used 
with the collision term associated to the grain boundaries. The background scattering was taken into account via the 
bulk relaxation time. The grains were considered to have columnar shape to be laterally bounded by planes. The 
grain boundaries were represented by two uncorrelated sets of (5-potentials along each lateral direction. The following 
expression was found for the conductivity of films 



(T (To 1 ~ 3q;o 



Qfo + In (1 + 1/ao) 



(30) 



where ctq is the conductivity of single crystalline material and the parameter ag is defined as 



"0 = 7; — -A 15- 31) 

DuJt- 1 — 7c 

Here Wi? is the Fermi velocity in the homogeneous material, D is the mean grain size, and TZ is the reflection coefficient 
of electrons on the grain boundary. It is important to stress that a depends only on the mean grain size; dependence 
on the standard deviation from the mean value can be neglected 41, 75j. Variation of the dc conductivity one can 
connect with the effective relaxation frequency, uj^^^ , by the relation 



Comparing it with Eq. H30I) one can find uo'^^^ . However, this simple concept does not work at frequencies where the 
conductivity demonstrates the dispersion. 

Recently the same approach was generalized to the case of optical conductivity cr (w) |4ll |. which is converging to 
the dc conductivity in the low- frequency limit, and the empirical parameter TZ was extracted from comparison of the 
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model with the measured reflectivity of gold films of different grain size. To find the optical conductivity the authors 
|4l| solved a linearized Boltzmann transport equation taking into account the anomalous skin effect as well. So, in 
general, the conductivity is nonlocal. It means that it depends not only on frequency but also on the wave vector. It 
was found that the anomalous skin effect gives only very small correction to the reflectance of the films. To simplify 
the model we will neglect the nonlocal effect. Recently the anomalous skin effect was discussed in connection with 
the Casimir force [l^l- It was found that the effect gives only a minor contribution to the force and for this reason it 
could be estimated separately for a homogeneous material neglecting the grain boundaries. 

In the local limit for the optical conductivity the expression similar to Eq. l|Sn|l was found [ill: 



a (uj) — (Jd {^) 1 — 3q; 



- — a + a'^ln(l + l /a) 



(33) 



where the conductivity in the Drude model ud (w) and the function a — a {u) are defined as 



The longitudinal component of the electric field, which is perpendicular to the film surface, does not feel the grain 
boundaries because the grains are columnar. Therefore, Eq. (|33|l is true only for the transverse field component. For 
this reason one has to distinguish two dielectric functions in the model: the longitudinal function ei {lj) = e (w), which 
coincides with the dielectric function of single crystalline gold, and the transverse function et {to). The latter one can 
be presented in the following form |4l| |: 



et (w) = e (tj) + 3 



a + a^ ln(l + l/a) 



(35) 



Using these dielectric functions Sotelo et al. [4l| predicted the film reflectivity that was compared with the measured 
reflectivity of the films with different mean grain size D. The experimental data were fitted in the visible and near 
infrared range to find the only parameter TZ. This parameter was estimated as 7?. = 0.650 ± 0.025. The results for 
different grain sizes all fall within the indicated errors. Prediction of the model was compared with the reflectivity in 
the infrared and excellent agreement was found. 



B. The Casimir force as a function of the voids fraction and the grain size 

To estimate influence of the film preparation procedure on the Casimir force, we use two parameter model taking 
into account the main reasons for the deviation of the optical properties from the bulk material. Voids in the film 
are described with the effective medium approximation (see Eq. H19|) l where the volume fraction of voids, /y, is 
considered as a free parameter. The other source of significant variation of the optical properties is the electron 
scattering on the grain boundaries. As was shown above this effect is completely defined by the mean grain size D 
if the reflection coefficient of electron, TZ, is fixed. The grain size D is the second free parameter which can be easy 
measured with AFM for any given film. Because the contribution of each effect, voids or grains, to the Casimir force 
is estimated on the level of a few percents, they can be considered separately. To make the calculations we have to 
fix the dielectric function of single crystalline gold. Instead of linear interpolation between e" (wi) and e" [ujq) in the 
Drude tail region, as was done for the upper limit evaluation, more realistic quadratic interpolation on the linear plot 
is used to fit the minimum e" (w) that happens at w « 2 eV for all the data. The real part of dielectric function, e' (w), 
does not show any special features in the Drude tail region. That is why for e' (w) the Drude function is used below 
(jJi and the handbook data 0,1^ above uji. Defined in this way function e' (u) shows much smaller discontinuity 
at w = wi than in the handbook between Dold and Mecke data [32'| and Theye data jsj . It should be stressed that 
these details in the dielectric function behavior can infiuence the force only on the level of 0.1%. 

It is relatively easy to calculate the change in the force due to variation of the void fraction in the film. It is 
assumed that the void fraction is small and one can use Eq. (|20|l for the effective dielectric function (e(cj)). This 
function is substituted in the dispersion relation Q to find {e(i()). The effective dielectric function at imaginary 
frequencies is used then in the Lifshitz formula to find the Casimir force. The relative correction to the force, AF/F, 
as a function of distance between bodies for fy = 0.05 is shown in Fig. |Slfor the sphere-plate geometry (dashed line) 
and plate-plate geometry (solid line). It was found that the dependence on the void fraction can be considered as 
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FIG. 8: The relative correction to the force due to presence of voids in the film as a function of separation. The results are 
given for the void fraction fv = 0.05. The dependence of the correction on fv is approximately linear. 

linear at least up to fv = 0.15. Hence, the graph can be used to estimate the correction for any reasonable value of 
fv- The correction is given for the case when the temperature effect is negligible (plasma model prescription for the 
n = term). The thermal correction can be estimated using Table Hll 

As was discussed in Sec. IIII Al all unannealed samples demonstrate presence of voids. The void fraction for the 
best unannealed sample evaporated with e-beam on NaCl substrate was estimated as fv = 0.04 assuming that 
annealed Theye films had no voids. For Johnson and Christy films 143 fv was estimated as 0.09 and the same value 
is realized for e-beam evaporated films in Wang et al. experiment [43 ■ We expect that fv = 0.04 is the smallest 
void fraction in unannealed samples because this value is realized for the films growing epitaxially '57*1. In the AFM 
experiment at the minimal separation a = 63 nm the correction to the force for the best unannealed sample is 
estimated then as —1.8 %. In the MEMS experiment 0|, where the smallest separation was a = 260 nm, for the 
correction it is found —1.0 %. These values practically coincide with that presented in Ref. where the calculation 
was performed using as the perfect material the handbook data extrapolated to low frequencies according to Ref. 
[25j | . This coincidence shows once more that the force is much more sensitive to the variation of the Drude parameters 
than to any other details of the dielectric function. Here and in Ref. 0| the same value of the plasma frequency, 
Wp « 9 eV, was used for the perfect material. 

It is not so straightforward to find the correction to the force due to variation of the grain size. The reason of 
complication is that the effect is described by two different dielectric functions e; (w) and eAuj). In this case the 
reflection coefficients and Vp have to be expressed via the surface impedances Zs and ^ 1791 . Discussion of the 
impedance approach to the Casimir force evaluation rose controversy in the literature 30, 72. ^ 731 l80l Isil Is^ "s?] . The 
problem was connected with the use of the approximate Leontovich impedance ^80,, ^3J for the force calculation. 
This approach was criticized |8lj| on the basis that, in general, one has to use two different impedances one for each 
polarization and it was demonstrated that the right local impedances reproduced the Lifshitz formula. Recently the 
problem was settled finally in Ref. where the classical theory of nonlocal impedances .JTOJ was applied to the 
Casimir problem. It was demonstrated that the Leontovich impedance is a good approximation for the propagating 
fields but it fails to describe the evanescent fluctuations that give significant contribution to the force. Calculating 
the correction to the force we follow the method developed in Ref. 73] . 

The Casimir force can be calculated using the same Lifshitz formulas Q and ^ where the reflection coefffcients 
are expressed via the impedances: 



In general, the impedances are functions of frequency C and wave vector component along the plate q. They can be 
written via the dielectric functions |79j and analytically continued to the imaginary frequencies |73| : 



c-z,.jeTq^ 

C + z,^C+¥^ 



(36) 
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2C /• dk, 

TTcJ (CVc2)£t(*C) + fc'' 




(37) 



2C 



(38) 



where fc = y'K^+q^ is the absolute value of the wave vector. Because in our case the dielectric functions are local 
(does not depend on k), the integrals here can be calculated explicitly. It gives 



Zs = 



(39) 



^Jet + icq/Cf 



(40) 



The dielectric functions that enter Zs and Zp can be written at imaginary frequencies with the help of Eq. 1)3511 as 



£/ «) = e «) 



£t «) = £ «) - 3- 



C(C + ^r) 

where e (i^) is the dielectric function of perfect material and 



a + a^ ln(l + 1/a) 



(41) 



7^ 



(42) 



To calculate the dependence of the Casimir force on the grain size, the impedances (|39|) and (|40|) were substituted 
in the reflection coefficients (I36|l which were used then in the Lifshitz formulas Correction to the force as a 

function of the grain size is shown in Fig. ^ The results are given for the closest separations a — 63 nm in the AFM 
experiment 9] (dashed line) and a = 260 nm in the last MEMS experiment (solid line). The correction found 
here is somewhat larger than in Ref. [23|. For the grain size Z) = 50 nm it is estimated as -0.9% for the AFM and 
-1.2% for the MEMS experiments in contrast with the corresponding values -0.5% and -0.9% reported previously [28l| . 
The difference is explained by the larger value of the relaxation frequency for the perfect material, uor — 0.035 eV ^ 
used in Ref. [2|. 

It was already stressed that, in general, the effect of finite grain size is not reduced to the effective relaxation 
frequency which follows from Eqs. (|Sn|l . 



1 — -ao -I- 3Q;g — 3q;q In (l + ctg ^) 



(43) 



However, the calculation showed that the correction to the force with a good precision is given by the relation 



— ^Aia)"^^^^. (44) 

t UJr 

This relation holds true at least in the investigated region 20 < 13 < 200 nm and 60 < a < 300 nm. The function 
A{a) is shown in Fig. ^|for the sphere-plate geometry (dashed line) and for the plate-plate geometry (solid line). 
Therefore, the concept of the effective relaxation frequency is happened to be good in practice. 
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FIG. 9: The relative correction to the force due to electron scattering on the grain boundaries as a function of the mean grain 
size. The results are given for the smallest separations in the experiments g| (dashed line) and [T3 | (solid line). 
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FIG. 10: The correction due to electron scattering on the grain boundaries as a function of separation. The function A(a) is 
connected with the correction by the relation 14411 . 



The grain size is difficult to predict because it depends on many details of the deposition procedure. Definitely one 
can say only that for evaporated films the grain size increases with the film thickness h as In general, 

it is not true that the grain size is close to the film thickness. For example, the grain size for thermally evaporated 
film 'ssj with the thickness of 205 nm can be estimated from the presented image as 13 w 50 nm. The image of 120 
nm thick evaporated film |67| gives D « 40 nm. 

The discussion above shows that instead of the physical parameters defining the film optical properties such as the 
volume fraction of voids /y and the grain size D one can use the effective Drude parameters connected with /y and 
D by Eqs. (|22|) and (|43|l . These effective parameters can be extracted directly from the optical data in the infrared 
as was described in Sec. 1111 CI The Drude parameters are the most important information about the film that one 
has to know to make high precision prediction of the Casimir force. The details of the film dielectric response in the 
near infrared and visible range can change the force only on the level of 0.5% or smaller. 

Voids and grain boundaries are the main but not the only factors affecting the electron system in the metallic film. 
The effective mass of electron, m*, can be modified by the strains which exist in the film and change slightly the 
crystallographic parameters. It would modify the Fermi surface but the changes are not sufficient to account for the 
observed variation of Wp [3l|. Some minor contribution, however, is not excluded. Internal electron-scattering defects 
like impurities or local defects can affect the relaxation time. This mechanism becomes dominant over the grain 
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boundaries in alloys. For monoatomic material like Au one can hardly believe that the local defects (0-dimensional) 
can play more significant role than the grain boundaries (2-dimensional). The Drude tail is very sensitive to the 
details of the film preparation |^ |H3| . The mechanisms of this dependence are not well understood [sJl but we have 
seen that the Drude tail can change the force less that 0.5% at the smallest separation. 

VI. DISCUSSION 

A number of arguments were provided in this paper to show that the optical properties of deposited gold films 
differ from the properties of bulk material on the level which cannot be ignored in the precise evaluation of the 
Casimir force. The best properties demonstrate carefully prepared annealed films but all the experiments were made 
with unannealed films. The main physical reasons for the deviation from the bulk properties are the voids in the 
films reducing the density of conduction electrons and the electron scattering on the grain boundaries significantly 
increasing the effective relaxation frequency. It was demonstrated also that the Casimir force is much more sensitive 
to the material optical properties in the infrared and far infrared range than to the dielectric response in the near 
infrared and visible range. For this reason the precise values of the effective Drude parameters are the most important 
information that has to be extracted from the optical measurements for successful prediction of the force. 

A reliable way to get this information is the direct optical characterization of the films used for the Casimir force 
measurement. Ordinary ellipsometry in the visible range, 1.5 — 4.5 eV, can be helpful to see the quality of the 
film but cannot be used to get the Drude parameters. Both real and imaginary parts of the dielectric function in 
the wavelength range 2 — 30 fim can be routinely measured with the infrared spectroscopic ellipsometcrs (for recent 
publications using these apparatus see 86, 87, 88J). The data can be fitted with the Drude dielectric function to 
find the effective plasma frequency cUp and the relaxation frequency uj-r characterizing the film. Special attention has 
to be paid to the sample preparation procedure to reduce the Drude tail and roughness of the film. Up to date the 
best known method for metal deposition is the filtered arc deposition process. The gold films deposited with this 
method demonstrated the bulk-like behavior in the near infrared |44| without annealing. The method is superior to 
evaporation or sputtering due to significant fraction, 70-80%, of Au ions with the energy 40 — 50 eV in the plasma. 
For high precision experiments the sample preparation and characterization is as important as the measurement of 
the force itself. 

While there is no direct optical data for the films used in the Casimir force experiments, one can establish the upper 
limit on the force usin g th e dielectric function of single crystalline gold. This upper limit was discussed in Sec. IIVI 
It was already noticed l23| that the AFM experiment 9] gives at smallest separations the force which is larger than 
the upper limit. Much more detailed analysis of the dielectric function of perfect single crystalline gold carried out in 
this paper supports this conclusion. 

Let us be more precise. For illustration we consider only one experimental point of the closest approach in the 
AFM experiment, Jjut the same conclusions will be true for all the points with a < 67 nm. The first point presented 
in Fig. 5 of Ref. 9] corresponds a = 63 nm and F^xp = 492 pN. This and the other points of close approach were 
found from the original postscript figure presented in the preprint quant-ph/0005088 with the digitizing errors (0.5 
pixel) 0.15 nm in the separation and 0.4 pN in the force. One can compare F^xp with the upper limit on the force 
which was discussed in Sec. IIVI For the reduction factor it was found ri{63nm) = 0.445 that is equivalent to the 
force F(63nm) — 464 pN. The difference F (63 nm) — Foxp — —28 pN should be compared with the experimental 
error that was estimated as AFcxp — 3.5 pN j^]. 

Due to error in the absolute separation Aa = 1 nm |^ the first experimental point could correspond to a = 62 nm. 
In this case the upper limit is rj {62 nm) = 0.442 and F {62 nm) — 483 pN. It should be stressed that the theoretical 
uncertainty in the force because of Aa is [F (62 nm) — F (63 nm)] / F (63 nm) = 4.1%. It is considerably larger than 
the claimed 1% agreement between the theory and experiment This problem has been already emphasized before 
[2^. Even for a = 62 nm the upper limit is still 9 pN smaller than the measured force. 

In reality the theoretical force should be smaller than the upper limit due to deviation of the film optical properties 
from the properties of bulk material. The smallest volume fraction of voids was estimates in Sec. IV Bl as 4%. It 
reduces the absolute value of the force at a = 63 nm on 1.8%. The electron scattering on the grain boundaries 
reduces the force additionally. As the largest possible value of the mean grain size for the unannealed 87 nm-thick 
gold film we take D — 100 nm. According to Eq. H30() it corresponds to the resistivity 4.4 fiil ■ cm or equivalently 
to the effective relaxation frequency w^-^-'^ = 0.044 eV. The minimal reduction of the force due to finite grain size 
is estimated then as 0.5%. Therefore, in the very best situation the theoretical force is expected to be smaller than 
F (63 nm) (1 — 0.018 — 0.005) = 453 pN or 472 pN for a = 62 nm. The difference with the experimental value is 
striking. 

Recently the result of the AFM experiment |9| has been reanalyzed ||21|. Three of 30 scans were rejected due to 
excessive noise while the rest 27 scans were used to calculate the force and the standard deviation. Smaller value was 
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found for the standard deviation Ai^exp = 2.8 pN . The only point for the force was reported: at a = 62 nra the mean 
force was Ff.^p = 486 pN . Both position and the force are smaller than that in Ref. ^ but the authors did not make 
any comment in this respect. The problem with the uncertainty of absolute separation Aa was dismissed on the 
basis that the experiment is in a good agreement with the theory. The theoretical result for the force was found as 
in Ref. using the handbook data and the Drude parameters ujp — 9.0 eV , uJr = 0.035 eV to extrapolate to the 
low frequencies. We have seen in Sec. IIVI that this procedure gives the force very close to the upper limit. The films 
used in the experiments should have different Drude parameters. In this respect the authors |2l| stated that uip is 
determined by the properties of the elementary cell which cannot be influenced by the sample quality. This statement 
is wrong because, as shows Eq. (fTUjl . ujp depends on the concentration of the conduction electrons which, indeed, 
depends on the sample quality. The voids will reduce the electron concentration in the sample. Even more confusing is 
the discussion of the sample dependence of iOr . The authors introduce not imaginary but real characteristic frequency 
(^c/i — c/2a and assume that this frequency gives the main contribution to the force. As was discussed in detail in Sec. 
^low frequencies cu ^ ojch give the main contribution to the dispersion integral (O for e (iC). This fact was already 
stressed in a number of publications [2^ l2^ 12^ . l26l | . Because of this mistake the authors took into account some 
minor effects in the near infrared like electron-electron scattering that gave quadratic frequency dependence of the 
relaxation frequency [sij but missed very important and well established effect like significant increase of the film 
resistivity in comparison with the bulk material. 

The same analysis cannot be done for the latest MEMS experiment [3! mainly because of huge roughness correction. 
In contrast with the AFM experiment, where the roughness correction was smaller than 0.5%, in this experiment the 
roughness correction is as large as 25%. One can hardly believe that the calculations based on the local separation 
[T^ can be true for the shape dependent Casimir force with the precision on the level of 1%. In the experiment 0| 
one of the plates was covered with copper but in the present paper an extensive analysis of gold was made. For this 
reason one can say only that the sample effect should reduce the force on the level of 2%. It is much larger than the 
declared experimental errors. The roughness correction and sample effects do not allow us to support the conclusion 
[T^ on the possibility to distinguish between different approaches to the temperature correction. 



VII. CONCLUSIONS 

In this paper we analyzed the sample dependence of the Casimir force. It was already noted [27Ll8^ that the optical 
properties of metallic films used to measure the force can deviate significantly from the bulk material but here, for 
the first time, serious analysis of the effect based on the facts about gold films was undertaken. It was stressed that 
the behavior of the dielectric function of a metal in the far infrared, where direct optical data are inaccessible, is 
crucial for the precise evaluation of the force. The Casimir force depends on the dielectric function e connected 
with the observable function e" (oj) via the dispersion relation Q. Because e" (uj) is large at low frequencies, the far 
infrared range gives the most significant contribution in e and, therefore, in the force. For example, at ^ = 1 eV 
the contribution of the real frequency domain u; < 0.125 eV to e was estimated as 75%. 

Sample dependence of the optical properties of gold films was discussed in Sec. IIIII The data for 40 years of 
investigations have been collected and analyzed. These data clearly showed that significant difference between samples 
cannot be explained by the experimental errors and had to be attributed to genuine sample dependence of the optical 
properties. The mechanisms of this dependence have been clarified in the literature. Variation of e" (uj) in the visible 
range was connected with the volume fraction of voids fv in the films. The void fraction is especially significant in 
unannealed films. For the best unannealed films which start to grow cpitaxially on NaC'l substrate fy was estimated 
as 4% but in the worst cases it could be more than 20%. The fraction of voids is very important for the Casimir force 
because it defines the precise value of the plasma frequency. The concentration of conduction electrons in the sample 
decreases with the increase of fv and, as the result, ujp decreases. 

Analysis of the data available in the infrared range revealed the same tendency of the dependence on the sample 
preparation procedure. It was demonstrated that at w < 1 eV the data for real and imaginary parts of s (uj) could 
be reasonably well described with the Drude model. The fitting procedure allowed one to find the effective Drude 
parameters for the investigated samples (see Table Pi . Surprisingly, the handbook data in the range oj < 1 eV 
taken from the old paper by Dold and Mecke [s^l happen to be bad. The best fit gave ujp = 7.52 eV, lUt = 0.061 eV. 
It means that the films were poorly prepared and were full of voids and defects for the electron scattering. On the 
other hand bulk ylit(llO) and well annealed gold films demonstrated much better characteristics close to that typically 
used for the calculation of the Casimir force |25|: ojp = 9.0 eV, cUr = 0.035 eV. 

It was stressed that the force was measured between unannealed films for which one can hardly expect that the 
effective Drude parameters coincide with that for the bulk material. However, using the "ideal" dielectric function of 
gold one can establish the upper limit on the Casimir force. This dielectric function was carefully described in Sec. 
IIVI The upper limit happened to be very close to the force evaluated in Ref. . The use of nearly the same ujp and 
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ujr for "ideal" gold in this paper and in Ref. '25!| explains the coincidence of the forces. Thus, the result of Ref. 
should be considered rather as the upper limit on the force than as the force between real films. 

The expected deviations of the force from the case of perfect single crystalline gold were analyzed in Sec. An 
appreciable correction to the Casimir force was found due to the presence of voids in the film. Even for the smallest 
fraction of voids, /y — 0.04, the relative correction to the absolute value of the force in the AFM experiment 
was found to be -1.8% at the smallest separation a = 63 nm. For the MEMS experiment at a = 260 nm the 
correction was -1.0%. The other effect that was estimated was influence of the grain boundaries in the film on the 
electron scattering. This effect is the main reason for the increase of the film resistivity or the effective relaxation 
frequency. For the mean grain size of 50 nm (typical for 100 nm thick film) the corrections due to the finite grain size 
were estimated as -0.9% for the AFM experiment and -1.2% for the MEMS experiment at the smallest separations 
in each case. Both corrections cannot be ignored in the experiment where the force was measured with the precision 
on the level or better than 1%. A number of additional factors can affect the dielectric function. Not all of them are 
well understood but one can indicate strains in the film or some internal electron-scattering defects. All additional 
effects are expected to reduce the absolute value of the force less than 1%. 

The best way to get the precise values of the Drude parameters is the direct measurement of the dielectric functions 
of the film used in a Casimir force experiment. The most important frequency range where this function should be 
known is the infrared range lo < 1 eV. It can be done with the infrared spectroscopic ellipsometers. Special attention 
has to be payed to the film preparation procedure. 
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